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ABSTRACT 

Although there is considerable evidence supporting an ubiquitous magnetic field in solar/stellar pho- 
tospheres, its impact in the determination of abundances has never been quantified. In this work we 
investigate whether the magnetic field plays a measurable role for this kind of studies. To that end, we 
carry out simulations of spectral line formation in the presence of a magnetic field, and then use those 
profiles to derive the abundance of several atomic species (Fe, Si, C and O) neglecting the magnetic field. 
In this way, we find that the derived iron abundance can be significantly biased, with systematic errors 
up to 0.1 dex. In the case of silicon, carbon and oxygen their role is very marginal (errors smaller than 
0.02 dex). We also find that the effect of the magnetic field strongly depends on its inclination with 
respect to the observer. We show that fields that are aligned with the observer lead to an underestima- 
tion of the real abundance, whereas more inclined ones overestimate it. In the case of a mixture of fields 
with different inclinations these effects are likely to partly cancel each other out, making the role of the 
magnetic field even less important. Finally, we derive a simple model that can be used to determine the 
suitability of a spectral line when we wish to avoid the bias introduced by the neglect of the magnetic 
field. 

Subject headings: Sun: magnetic fields, abundances - Stars: magnetic fields, abundances 



1. INTRODUCTION 

From the very first determinations of solar and stellar 
abundances using ID semi-empirical atmospheric models 
(e.g. Lambert 1968, Lambert & Warner, B. 1968, Grevesse 
1968, Garz et al. 1969) to more recent values obtained 
from state-of-the-art 3D hydrodynamical simulations (e.g. 
Asplund et al. 2005a), the role of the magnetic field has 
rarely been considered. This would be strictly valid only 
if atomic transitions with zero Lande factors are used in 
the analysis. Unfortunately this has never been the case. 
The reason for this is that there are few magnetic insensi- 
tive spectral lines having accurate oscillator strengths. In 
the case of the Sun, the role of the magnetic field has been 
avoided by arguing that the FTS-disk center spectral atlas 
(Brault & Neckel 1987; Neckel 1999), which is the most 
common source to compare observed and simulated line 
profiles, was recorded around a quiet Sun region. 

The existence of significant magnetic flux in quiet Sun 
regions has passed unnoticed because this magnetic field 
organizes in patches of opposite polaritjQ over very small 
scales, leading to a cancellation of the polarization sig- 
nals. However, there is now strong evidence supporting 
the omnipresence of magnetic fields in regions previously 
thought to be void of them (Trujillo Bueno et al. 2004; 
Manso Sainz et al. 2004). The details about its actual 
strength and distribution are subject to debate, however. 

This is even more critical when disk integrated data is 
used in stellar abundances studies (e.g. Allende Prieto et 
al. 2002), in particular if the star is magnetically active 
(Ap-Bp types), as regions of strong magnetic field (i.e.: 
starspots) can have a large influence. 

Although the magnetic field will mainly affect the po- 
larization signals, it also has an impact on the intensity 



profiles (Stokes /), being this particularly important be- 
cause this effect on the intensity profiles adds up regard- 
less of its polarity. If not accounted for, this might lead 
to systematic errors in the abundance values derived from 
the fitting. The goal of this work is to asses the impor- 
tance of neglecting the role of magnetic fields. To that end 
we will perform simplistic simulations of intensity profiles 
of various important atomic elements in the presence of 
a magnetic field. In Section 2 we describe the synthesis 
code employed, the spectral lines used and we briefly re- 
view the Zeemann effect applied to Stokes /. In Section 3 
we study how the strength and inclination of the magnetic 
field taintes the inferred abundances of iron, silicon carbon 
and oxygen when the existence of this field is neglected. In 
Section 4 we derive a simple model that is able to quanti- 
tatively predict the commited error. Section 5 is devoted 
to studying if results from our simple modeling could differ 
significantly if the same investigation is performed using 
more realistic 3D MHD models. Section 6 summarizes our 
findings and anticipates possible future work. 

2. SPECTRAL LINES SYNTHESIS 
2.1. 



Synthesis code 

We have employed the SIR code (Ruiz Cobo & del Toro 
Iniesta 1992) to produce synthetic spectral lines. This code 
solves the radiative transfer equation in the presence of 
a magnetized plasma. Although SIR synthesizes the full 
Stokes vector, we will restrict ourselves to consider only 
the total intensity, Stokes /, as done in most abundance 
studies. In addition, SIR allows for the magnetic field vec- 
tor to be a function of the optical depth, but in this work 
we will consider it constant. 

The Harvard-Smithsonian Reference Atmosphere 

1 In spectropolarimetry, the term positive polarity refers to magnetic fields pointing towards the observer: inclination angles 7 < 90°, while 
negative polarities refer to magnetic fields pointing away from the observer: 7 > 90° 
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(HSRA; Gingerich et al. 1971) was used in our calcula- 
tions. Note that using this ID LTE semi-empirical model 
means that we study the impact of the magnetic field 
alone, that is, we assume that other important ingredients 
of the spectral line formation are already being accounted 
for. For instance, the magnetic field couples with the 
energy and momentum equation, resulting in a modifica- 
tion of the temperature stratification. In our analysis we 
are assuming that the correct temperature stratification 
and convective velocity fields are known (i.e.: obtained 
through realistic 3D simulations) and therefore we focus 
on the magnetic field. If this was not the case, small errors 
in the temperature or velocities would dominate over the 
neglect of the magnetic field. 

To account for the convective broadening of the spectral 
lines we use, unless otherwise specified, a macroturbulent 
velocity of 2 km s _1 and a microturbulent velocity of 1 km 
s _1 . We anticipate (see Section 5) that these values have 
no particular consequences in our discussion. 

2.2. Spectral line selection 

We have decided to focus our study on four important 
atomic elements: Fe I, Si I, C I and O I. Iron is particu- 
larly important because it shows a large number of atomic 
transitions at visible wavelengths, and therefore it is of- 
ten used to investigate solar and stellar atmospheres. In 
addition, it is commonly used to distinguish whether a 
star is first or second generation since heavier elements 
are only produced in Supernova explosions (Christlieb et 
al. 2002). Silicon is important in the solar context as it 
is used, together with iron, to compare with meteoritic 
abundances (Asplund 2000). Oxygen and carbon's impor- 
tance comes from being the third and fourth most abun- 
dant elements in the Universe, respectively. In addition, 
they have both been targeted as being responsible for the 
lowering of the solar metallicity (Asplund et al. 2005a), 
that has caused major discrepancies between helioseismic 
inversions and solar models (Castro et al. 2007). There- 
fore, it is worthwhile investigating whether magnetic fields 
could have something to say in this regard. 

In total we consider 57 spectral lines: 29 of Fe I, 15 of 
Si I, 9 of C I and 4 of O I. Their properties are summarized 
in Table 1. They have been adopted from Asplund (2000) 
and Asplund et al. (2000, 2004, 2005b). We have rejected 
those lines for which accurate collisional parameters (un- 
der the ABO theory; Anstee & O'Mara 1995; Barklem & 
O'Mara 1997; Barklem et al. 1998) were not found. Note 
that only one line, Si I 5665.555 A has a zero Lande fac- 
tor: g c g = 0. Several neutral iron and carbon lines are 
potentially very sensitive to magnetic fields: g e s ^ 2. 

2.3. Intensity profiles under the presence of a magnetic 

field 

Ignoring off-diagonal elements in the propagation ma- 
trix, the main contributor to the intensity profiles is r\\ 
(Wittmann 1974; del Toro Iniesta 2003): 

m = 1 + y sin 2 7 + X - [<f> b + cf> r ](l + cos 2 7 ) | (1) 

where 770 is related to the abundance of the element, the 
excitation potential of the lower level and the transition 



probability of the atomic transition (Landi DcgPInnocenti 
1976). 7 refers to the inclination of the magnetic field 
vector with respect to the observer. Since it appears as 
cos 2 7 and sin 2 7, its contribution is the same regardless 
of the polarity of the magnetic field (see Footnote 1). 
Treating the real Zeemann pattern as an effective triplet 
(J = 1 — > J = 0), the functions <p pi <p r and refer to 
the Voigt profiles for the n (AM = 0), blue (AM = -1) 
and red (AM = 1) components of the Zeeman pattern, 
respectively. The former is centered at Ao (central labo- 
ratory wavelength) , whereas the latter ones are shifted by 
an amount ±Xb with respect to the II component. 

A B = Cg cS B\ 2 (2) 

where C = 4.67 x 10~ 13 [A Gauss] -1 , B is the strength 
of the magnetic field (measured in Gauss), and g e g is the 
effective Lande factor of the spectral line. 

A number of remarks are in order. First of all, in the 
absence of a magnetic field (or if g g = 0), <p p = (j> r = <pb 
and therefore r\i = 1 + rjo<p p /2. Also, it is interesting to 
see that the blue and red components are always present 
independently of the orientation of the field. When the 
field is small enough the splitting is much smaller that the 
Doppler width of the spectral line, As <C A Ad, causing the 
spectral line to broaden. In the opposite case, > A Ad, 
they appear as two separate spectral lines with half the 
strength of the spectral line in the absence of magnetic 
field. Finally, the central II component vanishes for mag- 
netic fields aligned with the observer (7 = 0) causing the 
line core to desaturate. 



3. EFFECTS ON THE MAGNETIC FIELD IN THE 
ABUNDANCE OF FE, SI, C AND O 

We now proceed to calculate synthetic Stokes / pro- 
files of the spectral lines in Table 1. We use the numerical 
code and the atmospheric model described in Sect 2.1. We 
consider a magnetic field with a varying strength and in- 
clination: B = [0, 500] Gauss, 7 = [0, 90] deg. We consider 
also standard solar abundances: loge^ e = 7.45 (Asplund 
et al. 2000), loge s * = 7.51 (Asplund 2000), loge c = 8.39 
(Asplund et al. 2005b) and logeo = 8.66 (Asplund et al. 
2004). The resulting profiles are then fitted with differ- 
ent abundances but assuming that there is no magnetic 
field. For Fe I we use \oge F e = [7.30, 7.80] in steps of 0.01. 
For Si I and O I we consider log est — [7.48, 7.54] and 
logeo = [8.62,8.70], respectively, both in steps of 0.002. 
Finally, for C I we use logec = [8.32,8.46] in steps of 
0.005. 

From the comparison of the original profiles with mag- 
netic field and fixed abundance with those where the abun- 
dance varies and the magnetic field is neglected, we obtain 
a x 2 ~ log e curve for each spectral line. We therefore in- 
fer for each line an optimum abundance, as the one that 
minimizes the \ 2 curve. No other free parameters are con- 
sidered. An example of this process is presented in Figure 
1 using the spectral line Fe I 6136.994 A a magnetic field 
of B = 250 Gauss and 7 — 60°. The final inferred abun- 
dance is obtained as the mean of the best-fit abundances 
from all spectral lines of that atomic element. 
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B=250 [G]; 7=60 [deg]; Original log€ F ,=7.45 
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Fig. 1. — x 2 between the intensity profile of Fe I 6136.994 A 
with magnetic field (B = 250 Gauss; 7 = 60°) and an iron abun- 
dance of loge^? e = 7.45, and the intensity profile of the same spec- 
tral line without magnetic field, as a function of the abundance. 
The abundance that produces the best fit to the former profile is 
loge Fe ^ 7.49. 



Figure 2 presents the errors A log e = log e^t — log e rca i 
(as a function of the magnetic field strength and inclina- 
tion) introduced in the abundance of Fe I, Si I, C I and O I, 
when the magnetic field is not accounted for. As expected, 
the larger the field strength the larger the error. Neutral 
iron, with up to Alogep < 0.1, presents the largest devi- 
ations. For the rest of considered elements the magnetic 
field seems to have only a marginal effect: A log esi < 0.01 
and Aloge c ,o < 0.02.' 

We also find that vertical fields underestimate the cor- 
rect abundance, whereas the opposite happens for more 
inclined magnetic fields. This can be explained attend- 
ing to Sect 2.3 (Equation 1). For vertical magnetic fields, 
7 = 0°, the II component of the Zeemann pattern is ab- 
sent, leading to a desaturation of the core intensity and 
therefore requiring a smaller abundance to fit the line pro- 
file. An example with Fe I 5250.209 A is presented in 
Figure 3 (left panel). However, when the magnetic field 
is horizontal, 7 = 90°, the line mainly broadens and thus 
yieling a larger abundance. See example for Fe I 6200.313 
A in Figure 3 (right panel). 
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Fig. 2. — Errors in the inferred abundances Aloge = loge^t — 
logical when the role of the magnetic field is neglected, as a func- 
tion of the magnetic field strength for 4 different inclinations: 7 = 0° 
(dotted line; vertical magnetic field), 7 = 30° (dashed), 7 = 60° 
(dashed-dotted) and 7 = 90° (dashcd-triple dotted; horizontal mag- 
netic field). Solid lines indicate a mixture of field inclinations, in 
which profiles obtained with the previous 4 inclinations are aver- 
aged before inferring an abundance. 

We must also consider that different regions of solar and 
stellar atmospheres can posses magnetic fields with differ- 
ent inclinations. For instance, in the case of the solar gran- 
ulation, where the field is mostly horizontal in the granules 
(upflowing gas), but is generally vertical in intergranular 
lanes (downflowing gas). Since the inclination that mat- 
ters is with respect to the observer, this situation reverses 
as we move towards higher latitudes or towards the solar 
limb. Many other low and intermediate-mass stars also 
posses an outer convective layer. In addition, strong mag- 
netic concentrations also display a variety of inclinations: 
star-sunspots (umbra and penumbra), pores, network re- 
gions etc. Although they are normally avoided in the Sun, 
in the stellar case they certainly contribute to the observed 
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profiles. 

Since vertical and horizontal fields seems to have oppo- 
site effects in the derived abundance it is appropriate to 
study whether they can cancel each other out when a mix- 
ture of different inclinations is present. To study this effect 
we have carried out a similar experiment as the ones previ- 
ously presented, being now the difference that we average 
the intensity profiles obtained with 4 different inclinations, 
7 = 0, 30, 60, 90° (with the same field strength) before in- 
ferring an abundance for each spectral line. In this case, 
the errors in the retrieved abundance arc much smaller 
(see solid line in Figure 2), being only perceptible for the 
case of neutral iron: Alogepo < 0.02. This simulation 
is very simple in the sense that it does not include dif- 
ferent temperature stratifications, different velocity fields, 
etc. However, it helps to highlight that the aforementioned 
cancellation effect can indeed take place. This cancellation 
effect in Stokes / due to vertical and horinzontal magnetic 
fields, is similar in a way to the cancellation of circular 
polarization signals due to a mixture of magnetic fields 
pointing towards and away from the observer. 



case of Fe I only. It could be marginally important for 
Si I, C I and O I only if the magnetic field had a clearly 
preferred orientation. In stellar atmospheres, where we 
observe disk-integrated signals, this is unlikely the case. 
In addition, magnetic fields do not appear as a plausible 
source of error in the current controversy of the solar car- 
bon and oxygen abundance. 

4. A PHENOMENOLOGICAL MODEL 

The different behavior seen in Fe I as compared to the 
other three elements considered cannot be understood in 
terms of the magnetic sensitivity of the lines used, as 
most of them have similar Landc factors ranging from 
g e g = 1 — 2. The source of these differences must be there- 
fore thermodynamic. In this section we will develop a tool 
to differentiate whether a given spectral line of a particular 
atomic element is prone to yield unreliable abundances. To 
that end we assume that the error in the inferred abun- 
dances, that appears as a consequence of neglecting the 
magnetic field, is directly proportional to the changes in 
the line profile in the presence of a magnetic field: 



Fe I 5250.209 A 
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Fig. 3. — Filled circles are intensity profiles obtained with the 
nominal iron abundance of log ep e = 7.45, a magnetic field strength 
of B = 500 Gauss and inclinations of 7 = 0° (left panel) and 7 = 90° 
(right panel). Those profiles are then fitted without a magnetic field. 
When the original profile is affected by a vertical magnetic field, 
the best fit (solid lines) is obtained by underestimating the original 
abundance. The opposite occurs for horizontal magnetic fields. 

From these results we conclude that the effect of neglect- 
ing the role of the magnetic field in abundance determina- 
tions strongly depends on the atomic specie and spectral 
line under consideration, as well as the strength and incli- 
nation of the magnetic field. From the four atomic species 
considered, it appears to have a measurable effect in the 
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where we use maximum of the absolute value since the 
derivative changes sign as a function of wavelength. At- 
tending Equation 1, we will model the intensity profile as 
a combination of three Gaussian, where two of them are 
shifted by an amount ±\b with respect to the central lab- 
oratory wavelength Ao- The third Gaussian is centered at 
Ao and possesses twice the strength of the other two. 
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where Iq and AA refer to the core intensity and the 
Doppler width of the spectral line in the absence of mag- 
netic field. That is, 



lim /(A) = 1 - (1 - 7 )e" 

B^0 



(A-A )^ 



(5) 



Note that AA is related to the HWHM of the_ Gaus- 
sian profile by a constant factor: HWHM=A/log 2AA ~ 
0.832AA. Table 1 presents the values of Jo and AA ob- 
tained from a Gaussian fit to the intensity profiles in the 
absence of magnetic field for each spectral line. Finally, 
Equations (2)-(4) can be combined to write: 
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(6) 



This procedure is very similar to the one used in Cabrera 
Solana et al. (2005) but using Stokes I instead of the cir- 
cular polarization, Stokes V. Equation 6 allows to predict 
the effect of ignoring the magnetic field using: properties 
of the spectral line (regardless of the atomic specie) in the 
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absence of a magnetic field (I , AA, Ao), Lande factor g e g 
(given by the electronic configurations) and the strength of 
the magnetic field (used to calculate As). Note that the in- 
clination of the magnetic field does not play any role. The 
factor [cos 2 7 + 1] does not help to explain the fact that for 
vertical fields Aloge < 0, but Aloge > for horizontal 
fields. Therefore it is more appropriate to write: 



+a b k 



A log e = .4(7) max 

(A-A +A B )" 



1-/0 



2AA 2 



CgoffAg[(A- A + 



, . (A-A -A B )^ 

-(A-Ao-ABje 



(7) 



where the calibration constant A(j) can be determined 
for different field inclinations using our results in Section 
3. This is done in Figure 4, where we plot the left-hand 
side term of Equation 6 (obtained from Section 3) versus 
the right hand-side of the same equation (evaluated us- 
ing Table 1) for two limiting inclinations: 7 = 0° (upper 
panel) and 7 = 90° (lower panel) for all spectral lines. As 
it can be seen the correlation is good enough to empirically 
justify in first approximation, our assumption in Equation 
3. That is, the more a spectral line is affected by the mag- 
netic field, the more unreliable the derived abundance will 
be if the magnetic field is not considered. In addition, all 
spectral lines seem to follow a linear relation regardless of 
the atomic specie. The present model explains why Fe I 
is affected by the magnetic field more than Si I, C I and 
O I. As already mentioned, this is not due to the different 
Lande factors, but rather due to their different sensitivity 
to thermodynamic parameters: Iq and AA. Narrower and 
deeper spectral lines (AA and Iq small) such as the ones 
from Fe I, are more sensitive to the magnetic field than 
weak and broad spectral lines. 

It is also interesting to notice that the selected C I lines 
are broader and weaker than those of Si I and O I. Despite 
this, carbon is affected by the magnetic field as much as 
the other two (see Fig. 2). The reason is that the employed 
C I lines are located at larger wavelengths, Xq > Xq 1 '° , 
and therefore the Zeeman splitting is larger (Eq. 2 and 7). 

Since we use the same temperature stratification, the 
different Jo's and AA's among the selected spectral lines 
can only be due to differences in the excitation potential 
of the lower level Xh transition probability log gf, and col- 
lisional broadening parameters a and a. In first approxi- 
mation, the former two would be related to I whereas the 
latter two determine AA. 

The utility of the procedure described here lies in the 
fact that whenever another spectral line is considered (re- 
gardless of the atomic element), we can use this method 
to evaluate Equation 7, and therefore calculate its approx- 
imate position in Figure 4. This will give us an idea of the 
error introduced by ignoring the effects of the magnetic 
field, even if a numerical code that solves the radiative 
transfer equation in the presence of a magnetic field is not 
available. 



0.1 



0.05 



B = 500 [Gauss] y= [deg] 



? 0.00 



-0.05 



-0.1 





x Fe I 




X 


□ Si I 


- 


X 


a I 




XX 

XVm\ x xx 

&m. ~~ xx. x 

x xxS ~ ^x 


• c I 




. ^( 7 =0°) = -785.6 


5 ** - 


X 



0.00 0.25 0.50 0.75 1.00 1.25 

maxllei/3Bll [x10" 4 ] 



< 



0.4 
0.3 

0.2 - 

0.1 - 



B = 500 [Gauss] 7= 90 [deg] 



0.0 



^(7=90°) = 2897.6 



x* 

X - 

*x 



Fe I 
Si I 



□ 

A I 

r i 



-0.1 I 

0.00 



0.25 0.50 0.75 1.00 

maxllai/3Bll [xlO" 4 ] 



1 .25 



Fig. 4. — Error in the determination of element abundances when 
using individual spectral lines versus the derivative of its intensity 
profile with respect to the magnetic field. Red crosses show the 29 
Fe I lines from Table 1. Green squares indicate Si I. O I is denoted 
by blue triangles. Yellow circles show C I lines. A linear fit to the 
data points is indicated by the black dashed line. The slope of the 
fitted line is indicated by .4(7). Top panel: vertical magnetic field. 
Bottom panel: horizontal magnetic field. 

5. CONVECTIVE VS MAGNETIC BROADENING 

In our simple ID model the convective broadening is 
introduced through the macro and microturbulent veloci- 
ties, whereas in a more realistic 3D modeling this broad- 
ening naturally occurs when adding the profiles emerging 
from different regions (i.e.: granules and intergranules) in 
the solar atmosphere. At first glance it seems plausible 
for this convective broadening to mask the magnetic one, 
thus making the contribution of the magnetic fields even 
more negligible. However, we must take into account that 
the individual intensity profiles emerging from different re- 
gions of the solar atmospheres are already affected by the 
magnetic field. Therefore, when all intensity profiles are 
added up to produce the corresponding convective broad- 
ening, the fingerprints of magnetic field will still be visible. 

To prove this statement we have carried out a simple 
simulation of the process previously described. To that 
end we synthesize intensity profiles of Fe I 5247.05 A 
(geff=2) without magnetic field and no net line-of-sight 
velocity, but using a macroturbulent velocity of 2 km s _1 . 
The resulting profile is plotted in Figure 5 (I\ , thick dashed 
line). We then repeat this synthesis but now including a 
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vertical magnetic field, 7 = 0°, with a strength of 1000 
Gauss. The resulting profile is plotted in Figure 1 (I 2 , 
thick solid line). 

Next we perform a similar synthesis, where: a) we set 
the macroturbulent velocity to zero; b) we include a ver- 
tical field of 1000 Gauss and, c) we use a net line-of-sight 
velocity that changes from —4 km s _1 to 4 km s _1 . Each 
of the emerging profiles is indicated by one of the thin 
dashed lines in Fig. 1 (labeled as Ik)- With this we try to 
simulate the effect of having different structures affected 
by convective velocity fields. We then add all those pro- 
files using a weighting function with a Gaussian shape in 
order to mimic the macroturbulence. The resulting profile 
is indicated in Fig. 1 as J 3 (squares). 

*W / 7T -4 log 2 , n s 

where the FWHM of the Gaussian weighting function 
is equivalent to a macroturbulent velocity of 2 km s _1 . 
As it can be seen, ^3 = I 2 . Therefore, the broadening 
due to convective motions does not hide the broadening 
produced by magnetic fields. Our simulations do not con- 
sider that this magnetic field is probably different in gran- 
ules (upflows) and intergranules (downflows). We also use 
the same temperature stratification in all cases. Consid- 
ering a more realistic case (using results from 3D MHD 
simulations) might yield slightly different results, but the 
same basic idea will apply also in that case. In this exam- 
ple we have used a relatively strong magnetic field (1000 
Gauss) to facilitate a visual comparisons in Fig. 5, but it 
is worth mentioning that repeating the experiment with 
smaller fields leads to the same result. 



Fe I 5247.050 A 
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Fig. 5. — Simulated intensity profiles of Fe I 5257.050 AThick 
dashed line, I\ was obtained using a macroturbulent velocity of 
iinuc = 2 km s~ . In addition to this, we include a magnetic field 
with: B = 1000 Gauss and 7 = 0°, thus obtaining I2 (thick solid). 
Individual thin dashed lines, ij,, were obtained with ^ mac = 0, 
B = 1000 Gauss and 7 = 0° and a varying v\ oa from —4 km s — 1 to 
4 km s to simulate a convective velocity field. All ij. are added 
up into I4 according to Eq. 8 and 9. 

To provide further support to our argument we have 
recalculated the effects of the magnetic fields in the abun- 
dance (as done in section 3) but using a smaller macro- 
turbulence: v mac = 1 km s _1 . The resulting plots are 



identical to those in Fig. 2. This result might seem at 
odds with Eq. 7, since the new macroturbulent velocity 
changes the values of AA and Iq- However, if we repeat 
the calculation of the calibration curves in Fig. 4 no dif- 
ferences are observed in A log e as compared to the case 
with f mac = 2 km s _1 . In order for this to happen, the 
only possibility is that the horizontal axis changes. This 
results in a different calibration constant .4(7) but equal 
errors in the abundance. 

6. CONCLUSIONS 

Ignoring the broadening caused by the magnetic field 
when fitting the intensity profiles of spectral lines, may 
lead to an erroneous determination of atomic abundances 
in the Sun and other magnetically active stars. Although 
there have been previous works where the magnetic field 
has been considered (Kochukhov et al. 2004; Socas- 
Navarro & Norton 2007), to our knowledge this is the 
first systematic study aiming at quantifying the role of 
the magnetic field. Our results indicate that Fe I lines are 
more affected than lines from elements like Si I, C I or 
O I. We have shown that vertical magnetic fields lead to 
an underestimation of the real abundance, while horizon- 
tal fields tend to overestimate it. In a more real situation, 
where the spectral lines would receive contributions from 
magnetic fields with different inclinations, these two effect 
are likely to cancel each other, making the contribution 
of the magnetic field almost negligible with the exemption 
of perhaps Fe I. It is important to mention that the de- 
generacy between magnetic fields and abundances occurs 
only for small fields, since for very large magnetic fields, 
the spectral line is fully split in its different Zeemann com- 
ponents (see for example Nesvacil et al. 2004). We have 
also developed a phcnomcnological model that can used 
to determine if a particular spectral line is suitable for 
abundances studies that do not consider the effect of the 
magnetic field. 

Our analysis ID LTE analysis neglects the effect of hav- 
ing variations in temperature stratifications and/or con- 
vective velocities fields. Therefore, our results apply only 
if all those others possible sources of errors can be elimi- 
nated. If that was not the case, the error introduced by 
the magnetic field would be of second importance. This 
situation is highlighted by the work of Socas-Navarro & 
Norton (2007), who in spite of consistently considering 
the magnetic field, obtained strong discrepancies between 
the inferred abundances in quiet solar regions (e.g.: gran- 
ulation) and magnetic regions (e.g.: pores). The source 
of those discrepancies is therefore to be ascribed to small 
(< 100 K) errors in the temperature stratification. 

It would be ideal to repeat this work using realistic 3D 
MHD simulations (using several initial magnetic fluxes) 
in order to model more realistically the different temper- 
atures, velocities, field strength and inclinations present 
in solar and stellar photospheres. However, our work does 
provide a first hint that magnetic fields arc an unlike source 
of large errors in abundance determinations, unless a very 
particular spectral line or a very particular magnetic con- 
figuration is present. 

I wish to thank Eberhard Wiehr for bringing up the is- 
sue, during a talk I gave at the University of Gottingcn 
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in 2003, of using only spectral lines with zero Lande fac- 
tor in abundance determinations. That warning led to 
consider the effects of the magnetic field whenever those 
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were not available. Thanks also to Carlos Allendc Prieto 
and an anonymous referee for important suggestions and 
comments. 
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Table 1 

Summary of employed spectral lines; temperature parameter a and cross section for collisions with neutral 

HYDROGEN (iN UNITS OF BOHR'S RADIUS do) FROM ABO THEORY (BARKLEM ET AL. 2000); LANDE FACTOR G e ff IS 
CALCULATED UNDER LS APPROXIMATION FROM THE ELECTRONIC CONFIGURATIONS OF THE UPPER AND LOWER LEVELS. 



Atom 


\ r a i 




lo E(9J) 


a 


a/a 


9eff 


lo 


A A [mAJ 


Lower Level 


Upper Level 


T~l T 

le I 


I o o a n /i r 

4389.245 


0.052 


—4.583 


0.249 


217 


1.50 


0.338 


52.30 


5D3 


7F2 


be I 


r o a t nrri 

5247.050 


0.087 


—4.946 


0.253 


206 


2.00 


0.464 


59.36 


5D2 


7D3 


TTl t 

be I 


5250.209 


0.121 


—4.938 


0.253 


207 


3.00 


0.471 


59.21 


5D0 


7D1 


be I 


C A 1 O TAO 

5412.798 


/I AO A 

4.434 


1 TA 1 

— 1.761 


A OO A 

0.280 


AT1 

971 


A AT 

0.97 


A O A C 

0.845 


A A A A 

60.04 


■ ) f 1 A 

3(j4 


P-TT A 

5114 


be I 


rrnt C A A 

55z5.544 


/i OO A 

4.230 


1 AO A 

— 1.084 


A OOO 

0.238 


748 


1 C A 

1.50 


A CCA 

0.559 


A A AC 

66.65 


rrvA 

51)0 


5D1 


Fe I 


5701.544 


2.559 


—2.216 


0.237 


361 


1.12 


0.439 


70.38 


3F4 


3D3 


Fe I 


5784.658 


3.396 


—2.530 


0.244 


796 


1.87 


0.781 


64.04 


51 3 


5D4 


be I 


5956.694 


0.859 


—4.605 


0.252 


227 


0.70 


0.621 


64.91 


5F5 


7P4 


be I 


6082.710 


2.223 


—3.573 


0.271 


306 


2.00 


0.745 


65.48 


5P1 


3P1 


TTi„ T 

he I 


6136.994 


O 1 AO 

2.198 


O A C A 

—2.950 


A OA C 

0.265 


OO A 

280 


O A A 

2.00 


A C A C 

0.545 


TA OT 

70.27 


c no 

5F 2 


5D1 


TTl — T 

he I 


6151.618 


O 1 TA 

2.176 


O OA A 

—3.299 


A OA O 

0.263 


OTT 

277 


1 O O 

1.83 


A A OA 

0.639 


ATT A 

67.74 


c ~n 

5P3 


5D2 


be I 


/11 f7Q OOP" 

6173.335 


2.223 


—2.880 


0.266 


281 


2.50 


0.535 


71.16 


5P1 


5D0 


be I 


A AAA 10 

6200.313 


2.608 


—2.437 


0.235 


350 


1.50 


0.522 


73.08 


3F2 


3F3 


Fe I 


A O A 

6219.280 


2.198 


—2.433 


0.264 


278 


1.60 


0.436 


77.74 


5P2 


5D2 


Fe I 


A A A A O A O 

6240.646 


2.223 


—3.230 


0.272 


301 


1.00 


0.634 


68.97 


5P1 


3P2 


be I 


6265.133 


2.176 


—2.550 


0.261 


274 


1.58 


0.456 


76.98 


5P3 


5D3 


Fe I 


6271.278 


3.332 


—2.703 


0.247 


720 


1.50 


0.821 


69.31 


5F5 


7D5 


Fe I 


A O OA A 1 O 

6280.618 


0.859 


—4.387 


0.253 


223 


1.45 


0.565 


70.22 


5F5 


7F5 


Fe I 


6297.793 


2.223 


—2.740 


0.264 


278 


1.00 


0.507 


74.31 


5P1 


5D2 


be I 


A O O A A O C 

6322.685 


2.588 


—2.426 


0.238 


345 


1.50 


0.520 


75.02 


3F3 


3F4 


be I 


a /i 1 7A 

6481.870 


2.279 


—2.984 


0.243 


308 


1.50 


0.587 


73.46 


3P2 


5D2 


be I 


6498.939 


A A C O 

0.958 


A AAA 

—4.699 


A O C O 

0.253 


226 


1 OT 

1.37 


A A AO 

0.698 


A A OI 

69.81 


51 3 


TTPO 

71 3 


Fe I 


A C O 1 OI A 

6581.210 


1.485 


—4.680 


0.245 


254 


1.30 


0.859 


68.69 


3F4 


5F4 


Fe I 


/jrfiO O T A 

6593.870 


2.433 


—2.422 


0.247 


321 


1.15 


0.496 


79.91 


3H5 


r /~i c 

5C5 


Fe I 


a A AA 1 1 A 

6609.110 


2.559 


—2.692 


0.245 


335 


1.15 


0.589 


75.30 


3F4 


3C4 


be I 


6750.152 


2.424 


—2.621 


0.241 


335 


1.50 


0.542 


79.42 


3P1 


3P1 


TTU T 

be I 


a a a c o a c 

6945.205 


2.424 


O A OO 


A O A O 

0.243 


OOI 

331 


1 C A 

1.50 


A C 1 A 

0.519 


OA A 1 

84.01 


OT - ) 1 

3F1 


OTTO 


be I 


6978.851 


O AO A 

2.484 


O C A A 

—2.500 


A O A 1 

0.241 


O OT 

337 


1 C A 

1.50 


A C O A 

0.536 


O O AO 

83.43 


on a 

3P0 


OTT1 

3P1 


Fe I 


wwoo AAO 

7723.208 


2.279 


—3.617 


0.242 


304 


1.17 


0.794 


82.56 


OTOO 

3P2 


3D3 


Si I 


C A /I C A "1 O 

5645.613 


4.93 


—2.04 


0.223 


1791 


1.75 


0.748 


74.38 


3P1 


OOI 

3S1 


Si I 


5665.555 


4.92 


— 1.94 


0.222 


1772 


0.00 


0.708 


75.76 


3P0 


3P0 


Si I 


5684.484 


4.95 


— 1.55 


0.221 


1797 


1.25 


0.570 


82.41 


3P2 


OOI 

3S1 


C * T 

Si I 


C AAA /IOC 

5690.425 


A AO 

4.93 


1 TT 

— 1.77 


0.222 


1 TTO 

1772 


1 C A 

1.50 


A A A A 


TO A O 

78.42 


3P1 


OTT1 

3P1 


C * T 

Si I 


C TA1 1 A A 

5701.104 


A AO 

4.93 


1 AC 

— 1.95 


A OOO 

0.222 


1 TA T 

1767 


1 C A 

1.50 


A T1 A 

0.716 


TA A A 

76.06 


3P1 


OTO A 

3P0 


Si I 


C TA O /I A A 

5708.400 


4.95 


— 1.37 


0.222 


1787 


1.50 


0.509 


87.37 


3P2 


3P2 


Si I 


5772.146 


5.08 


— 1.65 


0.207 


2036 


1.00 


0.660 


81.16 


1 TO 1 

1P1 


1 OA 

ISO 


O * T 

Si I 


C TO A O O /I 

5780.384 


4.92 


—2.25 


0.228 


1691 


0.50 


0.819 


74.28 


3P0 


3D1 


Si I 


C TA O A ^7 

5793.073 


4.93 


— 1.96 


0.228 


1703 


1.00 


0.723 


77.01 


3P1 


3D2 


C * T 

Si I 


C TAT OCA 

5797.856 


A AC 

4.95 


1 AC 

— 1.95 


r\ 000 

0.223 


1 TC C 

1755 


1 1 A 

1.16 


A TOT 

0.727 


TT OT 

77.27 


3P2 


31)3 


C * T 

Si I 


C A /I O C A 1 

5948.541 


C AO 

5.08 


1 IO 

— 1.13 


A OOO 

0.222 


1 O A C 

1845 


1 AA 

1.00 


0.491 


A A C O 

96.52 


1P1 


11)2 


Si I 


7680.266 


5.86 


—0.59 


0.495 


2106 


1.00 


0.631 


126.0 


1P1 


1D2 


Si I 


TA IO O O A 

7918.384 


5.95 


—0.51 


0.232 


2933 


0.75 


0.643 


141.8 


3D1 


3F2 


Si I 


TA OO O A O 

7932.348 


5.96 


—0.37 


0.235 


2985 


1.00 


0.610 


150.3 


3D2 


3F3 


Si I 


TATA OAT 

7970.307 


5.96 


— 1.37 


0.232 


2927 


0.92 


0.883 


117.0 


3D2 


3F2 


C I 


T1 11 A O A 

7111.459 


O O A A 

8.640 


1 AT A 

— 1.074 


A OIO 

0.313 


1 O A O 

184z 


A TC 

0.75 


A A A 1 

(J. 941 


101 1 
121.1 


OlA1 


OTPO 

31 2 


C I 


T1 10 1 TA 

7113.179 


O O A T 

8.647 


A TAO 

—0.762 


A 1 A 

0.314 


1 O C T 

1857 


1 10 
1.12 


A AAO 

0.903 


IOC A 

125.4 


OT\0 

3D3 


OTT 1 A 

31 4 


C I 


9603.036 


7.480 


-0.895 


0.236 


561 


2.00 


0.752 


175.4 


3P0 


3S1 


CI 


10753.976 


7.488 


-1.598 


0.238 


532 


2.00 


0.875 


174.5 


3P2 


3D1 


CI 


11777.546 


8.643 


-0.490 


0.271 


746 


0.92 


0.854 


214.1 


3D2 


3F2 


CI 


12549.493 


8.847 


-0.545 


0.293 


863 


1.50 


0.880 


228.0 


3P0 


3P1 


CI 


12562.124 


8.848 


-0.504 


0.293 


863 


1.50 


0.874 


230.6 


3P1 


3P0 


CI 


12569.042 


8.848 


-0.586 


0.293 


863 


1.50 


0.885 


226.3 


3P1 


3P1 


CI 


12581.585 


8.848 


-0.509 


0.293 


862 


1.50 


0.874 


230.9 


3P1 


3P2 


01 


6158.176 


10.74 


-0.30 


0.322 


1915 


1.58 


0.976 


98.30 


5P3 


5D3 


01 


7771.944 


9.15 


0.37 


0.234 


452 


1.33 


0.734 


134.04 


5S2 


5P3 


OI 


7774.166 


9.15 


0.22 


0.234 


452 


1.92 


0.761 


129.76 


5S2 


5P2 


OI 


7775.388 


9.15 


0.00 


0.234 


452 


1.75 


0.800 


124.15 


5S2 


5P1 



